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We calculate the important next-to-leading-order contributions to the B — * ttK, nir decays from 
the vertex corrections, the quark loops, and the magnetic penguins in the perturbative QCD ap- 
proach. It is found that the latter two reduce the leading-order penguin amplitudes by about 10%, 
and modify only the B — > nK branching ratios. The main effect of the vertex corrections is to 
increase the small color-suppressed tree amplitude by a factor of 3, which then resolves the large 
difference between the direct CP asymmetries of the B° — > ir^K and B — > tt°K modes. The 
puzzle from the large B — > branching ratio still remains. 
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I. INTRODUCTION 

The B factories have accumulated enough events, which allow precision measurements of exclusive B meson 
decays. These measurements sharpen the discrepancies between experimental data and theoretical predictions 
within the standard model, such that some puzzles have appeared. The recently observed direct CP asymmetries 
and branching ratios of the B — > ttK, nn decays [J, 



A CP (B° - 
A C p(B ± - 
B(B° - 
B(B° 



TV K ) = 
• 7T + 7T~) = 



(-11.5 ± 1. 

(4 ± 4)% , 
(5.0 ±0.4) x 10~ 6 
7r 7r°) = (1.45 ±0.29) x 10 



-6 



(1) 



are prominent examples. The expected relations Acp(B° — ► ir^K^) Aqp{B ± — > 7r°ii' ± ) and B(B° — ► 
7r + 7r _ ) 3> B(B° — ► 7r°7r°) obviously contradict to the above data. In this work we shall investigate the 
indication of Eq. QJ, and study whether they can be accommodated in the perturbative QCD approach 0,0. 

To explain these puzzles, it is useful to adopt the topological-amplitude parameterization Q for these decays. 
The most general parameterization of the B — > 7T7t decay amplitudes is written as 



A(B Q - 
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(2) 



where T, C , P, and P ew stand for the color-allowed tree, color-suppressed tree, penguin, and electroweak 
penguin amplitudes, respectively, and (f>2 is the weak phase defined later. The counting rules in terms of powers 
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of the Wolfenstein parameter A ~ 0.22 are then assigned to various decay amplitudes 0,0,0- The amplitudes 
in Eq. obey the counting rules in the standard model, 



T T T 



A 2 



(3) 



Therefore, the B° — *■ 7r°7r° branching ratio is expected to be of 0(A 2 ) of the B° — > ir + ir~ one. However, Eq. (JIJ 
shows that the former is about of O(A) of the latter. 

The B — ► 7rif decay amplitudes are written, up to 0(A 2 ), as 



A(B+ - 




V2A(B+ - 




A(B° - 


> 77- K+) 


V2A(B° - 





P' 



P' 
r' 
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P 7 + P 7 



P' 



P' 

e-u 

p> 



p' 



(4) 



where the notations T' C , P' and P' ew bear the meanings the same as for the B — > irir decays but with 
primes, and the weak phase 4>3 is defined via the Cabibbo-Kobayashi-Maskawa (CKM) matrix element V u b = 
\V u b\ exp(—i(j> 3 ) @- These amplitudes obey the counting rules, 



pi ' pi ' pi 



A 



(5) 



The data Acp(B° — > it^K ) w — 11% implies a sizable relative strong phase between T 1 and P', which verifies 
our prediction made years ago using the PQCD approach 2j. Since both P' ew and C are subdominant, the 
approximate equality for the direct CP asymmetries Acp(P ± — > ir^K^) w Acp(P° — ■» ^J^) is expected, 
which is, however, in conflict with the data in Eq. |T|l dramatically. 

It is then natural to conjecture a large P' ew 0, [T(J, [U [13 , which signals a new physics effect, a large C 
[l3L IbH IT5I Il6| . or both 0,0] in the B — > 7ri^ decays. The large C" proposal seems to be favored by a recent 
analysis of the P — * ttK, tttt data based on the amplitude parameterization ^3|- Note that the current PQCD 
predictions for the two-body nonlcptonic B decays were derived from the leading-order (LO) formalism. While 
LO PQCD implies a negligible C, it is possible that this supposedly tiny amplitude may receive a significant 
subleading correction. Hence, before claiming a new physics signal, one should at least examine whether the 
next-to-leading-order (NLO) effects could enhance C sufficiently. In this paper we shall calculate the important 
NLO contributions to the B — > nK, tttt decays from the vertex corrections, the quark loops, and the magnetic 
penguins. The higher-power corrections have not yet been under good control, and will not be considered here. 
We find that the corrections from the quark loops and from the magnetic penguins, being about 10% of the LO 
penguin amplitude, decrease only the B — > ttK branching ratios. The vertex corrections tend to increase C 
by a factor of 3. This larger C leads to nearly vanishing Acp(P ± — > Tr ^*) without changing the branching 
ratios, which are governed by P' . The B — > irK puzzle is then resolved. 

The other NLO corrections, mainly to the B meson transition form factors, can be eliminated by choosing 
an appropriate renormalization scale /1 ~ y/Arrib, A being a hadronic scale and rru the b quark mass. This 
observation follows the well-known Brodsky-Lepage-Mackenzie (BLM) procedure [l!|, in which the scale fi is 
determined in the way that the vacuum polarization effects are absorbed into the coupling constant a s (/i). It 
has been demonstrated with this procedure that NLO corrections to many exclusive processes are minimized to 
some extent ^^|. Taking the simple pion form factor as an example, the BLM scale has been found to be of the 
order of the invariant mass of the hard exchanged gluon. The choice of \i proposed in the PQCD approach [2(| 
is basically in agreement with this procedure: the argument fi of the coupling constant is set to the invariant 
masses of internal particles, which are of 0(\/Amb) for the B meson transition form factors 0, |2^, Hj|- A 
general feature of the BLM scale is that it is always much lower than the external kinematic variable, implying 
that the smallness of the coupling constant is not the only condition for the applicability of perturbation theory. 

As mentioned before, the observed branching ratio B(B Q — > 7r°7r°) w 1.5 x 10~ 6 is much larger than the LO 
PQCD prediction ~ 10~ 7 0,|2J|. The prediction from QCD-improved factorization (QCDF) has the same order 
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of magnitude [25| . Since this mode involves a subdominant color-suppressed tree amplitude as shown in Eq. @ , 
a larger C certainly helps to resolve the B — > irir puzzle. We also compute the NLO corrections to these decays, 
and find the similar reduction from the quark loops and the magnetic penguins, which are about 10% of the 
LO penguin amplitude P. Since P is subdominant, the B° — * 7r =F 7r ± and B — » n^n branching ratios almost 
remain the same. The enhancement of C from the vertex corrections, leading to B(B° — > 7r°7r°) m 0.3 x 10~ 6 , 
is still not sufficient to account for the data. We point out that any new mechanism, introduced to resolve this 
puzzle, must survive the constraint from the tiny observed branching ratios [lj, 

B(B° A°A°) = x 10~ 6 , 

B{B° -> p°p°) < 1.1 x 10- 6 . (6) 

The leading PQCD predictions for B(B° -> K°K°) and for _BTB -> [HIH have been consistent with 

the experimental data. The proposals of final-state interaction 29] and of the charming penguin in soft-collinear 
effective theory (SCET) [3(j have not yet been applied to the B° — > p°p° decay. 

We review the LO PQCD predictions for the B — > irK, irir decays, including those for the mixing-induced 
CP asymmetries, in Sec. II. The vertex corrections, the quark loops, and the magnetic penguin amplitudes are 
computed in Sec. III. We perform the numerical study in Sec. IV, where the theoretical uncertainty is also 
analyzed. Section V is the conclusion. The explicit factorization formulas for the various topologies of decay 
amplitudes are collected in Appendix. 



II. LEADING-ORDER PREDICTIONS 



The effective Hamiltonian for the b — > s transition is given by [31j , 



10 



i=3 



(7) 



with the Fermi constant Gf = 1.16639 x 10 5 GeV 2 , and the CKM matrix elements V. The four-fermion 
operators are written as 
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(8) 



with the color indices i, j, and the notations (q'q')v±A = q'luO- ± 7s)?'- The index g' in the summation of the 
above operators runs through u, d, s, c, and b. The effective Hamiltonian for the b — > d transition is obtained 
by changing s into d in Eqs. Q and JSJ. 

According to Eq. Q, the amplitude for a B meson decay into the final state / through the b — ► s(d) transition 
has the general expression, 



A(B - /) = K* 6 K s ( d) + + TO s(d) Af 

For / = 7rA", the amplitudes A^, A^, and A„ K are decomposed at LO into 



(9) 



.4 



(u) 



rf - f K F e +M e + UFeK + M eK + f B F a + M a , 

A {c) - 

A^k = -{fi<F e p + M P + UF p K +M p K + f B F p +M p a ) 



(10) 
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where fs (/if, fir) is the B meson (kaon, pion) decay constant, F e (A4 e ) the color- allowed factorizable (non- 
factorizable) tree emission contribution, F e x {M-eK) the color-suppressed factorizable (nonfactorizable) tree 
emission contribution, F a (M a ) the factorizable (nonfactorizable) tree annihilation contribution, and those 
with the additional superscripts P the contributions from the penguin operators. For / = ttit, the amplitudes 
•A^m, -Attttj and Anl are decomposed at LO into 

A%> = f w F e + M e + f B F a + M a , 
A^ — 

•^irir — u ' 

A% = - (UF e p + M p e + f B F p + M p a ) , (11) 

where we do not distinguish the color-allowed and color-suppressed contributions. 

The factorization formulas for the various contributions to each B — ► irK, irir mode are collected in Tables UJ 
andHTI and in Appendix, whose dependence on the Wilson coefficients has been made explicit. We define the 
standard combinations, 



ai(fi) = C 2 (fx) + 
a 2 {p) = Ci{pb) + 



N c 

Nr. 



Oifr) = Ci(n) + , * = 3~10, (12) 

where the upper (lower) sign applies, when i is odd (even). The coefficients a and a' in Tables HI and ITU besides 
a\ and a 2 given above, are then written as 
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(13) 



With the amplitude _4(i? — * /) being computed using Eq. we derive the two-body nonlcptonic B meson 
decay rates and CP asymmetries. The former are given by 

T(B^f) = ^-\A{B^f)\\ (14) 

where ms is the B meson mass. The time-dependent CP asymmetry of the B° — > tt°Ks mode is defined as 

, , R o w o w v _ B(B°(t)^7r K s )-B(B (t)^K°Ks) 
A CP (B (t) - tt = ^o^^o^ + ^o^^) 

= ^ojcs cos(AM d i) + S 7r o Ks sin(AM d t) , (15) 

with the mass difference AMj of the two S-meson mass eigenstates, and the direct asymmetry and the mixing- 
induced asymmetry, 

. \K°K 3 ? - 1 Q _ 2Im(A 7r o Ks ) 

A ir°K s - TTTx 12 ' bir°K s - ~T~~~T\ [I > U b j 

respectively. The B° — ► ir°Ks decay has a CP-odd final state, and the corresponding factor, 

pi _ pi _ f]l p -i<t>3 
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where the weak phase 0i is defined via V t d — \V t d\ exp(— i<j)\). The time-dependent CP asymmetry of the 
— > 7T+7T mode is defined by 



A CP (B (t)^TT+7r-) 



B(B°{t) -» TT+7T-) - B(B°(t) -» TT+TT 



B(B°(t) -» 7T+7T-) + B(B°(t) -» 7T+7T- 

= Ajtt,. cos(AM d i) + sin(AM d i) , 
with the direct asymmetry and the mixing-induced asymmetry, 

_ 2Im(A 7r7r ) 



respectively, and the factor, 
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TABLE I: B — > 7rif decay amplitudes, whose factorization formulas are presented in Appendix. 
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where the weak phase 4>2 comes from the identity 4>2 = 180° — <j>\ — 4>3- In addition, the direct CP asymmetry 
for a charged B meson decay B + — ► f(B~ — > /) is dchncd by 

= B(B-^f)-B(B+^f) 
CP B(B~ -» /) + -»• /) ' 1 j 

The PQCD predictions for the branching ratios and the CP asymmetries of the B — > 7rif , 7T7t decays in the 
NDR scheme are listed in Tables Hill VII Using the LO and NLO Wilson coefficients, we obtain the values in the 
columns labelled by LO and LOnlowCj respectively. It is noticed that some of the NLO Wilson coefficients, 
like C5, diverge at a low scale. To derive the above tables, we have frozen the values Cj(//) at Cj(//o = 0.5) GeV, 
whenever fj, runs to below the scale /Xo, since the renormalization-group (RG) evolution is not reliable for \x < fi . 
Note that /iq is also the scale, which sets the starting point of the RG evolution of the Gegenbaucr coefficients in 
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TABLE II: £? — > 7T7t decay amplitudes, whose factorization formulas are presented in Appendix. 
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the meson distribution amplitudes [32j • We have kept the corrections in higher orders of the electroweak coupling 
a to the Wilson evolution, which were neglected in |3^|. Because the effect of the NLO Wilson coefficients is to 
enhance the penguin amplitude, the branching ratios of the penguin-dominated B — > irK modes increase, and 
the magnitudes of the direct CP asymmetries decrease a bit accordingly. As shown in Eq. J5J), The enhanced 
penguin amplitude P, being destructive to the color-allowed tree amplitude T, renders the B — > tttt branching 
ratios vary toward the direction favored by the data. The larger subdominant penguin amplitude also increases 
the magnitudes of the direct CP asymmetries in the B — » tttt decays due to the stronger interference with the 
dominant tree amplitudes. 

As stated before, the LO PQCD predictions for the B — > irK branching ratios are consistent with the data, 
viewing the range spanned by the columns LO and LOnlowc m Tabic 11111 However, the prediction for the 
B° — > 7T 7r° branching ratio is too small compared to the measured value. Those for the direct CP asymmetries 
of the B — > wK, tttt decays, except Acp{B ± — > -k^K^), are all in good agreement with the data as shown in 
Table ITVl The LO direct CP asymmetry of the B° — > 7r°7r° mode differs in sign from the result obtained in |3|], 
because we have employed the different pion distribution amplitudes (see Sec. IIVII . It simply implies that the 
theoretical uncertainty for the modes with tiny branching ratios is huge. Note that the predictions from QCDF 
plf for the direct CP asymmetries usually have signs opposite to those from PQCD. It has been realized that 
the set "S4" with non-universal parameters, such as the different annihilation phases for the B — ► PP, PV, 
and VP decays, must be adopted in order for QCDF to accommodate the data |33l l34l l35l l36| . The above two 
discrepancies associated with the B° — > 7r°7r° branching ratio and with the i? ± — > it^K^ direct CP asymmetry 
lead to the puzzles mentioned in Introduction. We prepare Table [V] for the various topological amplitudes, 
whose definitions are referred to 6]. The values in the columns LO and LOnlowc follow the power counting 
rules in Eqs. and JSJ exactly, explaining why the B — ► TtK, 7T7t puzzles appear. 





Mode 


Data [1] 


LO 


LOnlowc ±VC ±QL ±MP 


B ± 




24.1 ± 1.3 


17.0 


32.3 


31.0 34.2 24.1 


B ± 


iv K± 


12.1 ±0.8 


10.2 


18.4 


17.4 19.4 14.0 


B° 




18.9 ±0.7 


14.2 


27.7 


26.7 29.4 20.5 


B° 




11.5 ± 1.0 


5.7 


12.1 


11.8 12.8 8.7 



94 cj+13.6 (+12.9) 
8.1 (- 7.8) 

,o q +io.o(+ 7.0) 

5.6 (- 4.2) 
Of! n+15.6(+11.0) 
ZU - y - 8.3 (- 6.5) 
g j + 5.6 (+ 5.1) 



B° -nr^-ir* 5.0 ±0.4 7.0 6.8 6.6 6.9 6.7 6.5^^^^ 
B ± ^7r ± 7r° 5.5±0.6 3.5 4.1 4.0 4.1 4.1 4.0+ ^ 
1.45 ±0.29 0.12 0.27 0.37 0.29 0.21 0.29+^° 



TABLE III: Branching ratios in the NDR scheme (xl0~ 6 ). The label LOnlowc means the LO results with the NLO 
Wilson coefficients, and ±VC, ±QL, ±MP, and ±NLO mean the inclusions of the vertex corrections, of the quark loops, 
of the magnetic penguin, and of all the above NLO corrections, respectively. The errors in the parentheses were defined 
in the context. 





Mode 


Data [1] 


LO 


LOnlowc 


±vc 


±QL 


±MP 


±NLO 


B ± 




-0.02 ±0.04 


-0.01 


-0.01 


-0.01 


0.00 


-0.01 


0.00 ± 0.00 (±0.00) 


B ± 




0.04 ± 0.04 


-0.08 


-0.06 


-0.01 


-0.05 


-0.08 


n m +0.03 (+0.03) 
U ' U -0.05 (-0.05) 


B° 
B° 


. ^0 7^0 


-0.115 ±0.018 
0.02 ±0.13 


-0.12 
-0.02 


-0.08 
0.00 


-0.09 
-0.07 


-0.06 
0.00 


-0.10 
0.00 


n nq +0 06 (+ 04 ) 

u ' ua -0.08 (-0.06) 
n n? +0.03 (+0.01) 
U ' U ' -0.03 (-0.01) 


B° 


— > 7r :F Tr ± 


0.37 ±0.10 


0.14 


0.19 


0.21 


0.16 


0.20 


n 18 +0.20 (+0.07) 
U °-0.12 (-0.06) 


B ± 




0.01 ±0.06 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 ± 0.00 (±0.00) 


B° 


- rr'V 


n os+0-40 


-0.04 


-0.34 


0.65 


-0.41 


-0.43 


n fi o+0-35 (+0.09) 
U ' D,5 -0.34(-0.15) 



TABLE IV: Direct CP asymmetries in the NDR scheme. 



After obtaining the values of the various topological amplitudes, we compute the mixing-induced CP asym- 
metries through Eqs. (|T7I) and Since C is of 0(A 2 ) compared to P', it is expected that the LO PQCD 
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Topology LP LOnlowc +VC +QL +MP +NLO 

P' 36.6e* 2 ' 9 50.6e 129 49.6e i2 ' 9 52.1e i29 43.7e l28 44.1e l2 ' 9 
T" 6.9e l0 ° 6.6e' 6.6e l01 6.6e l0 ° 6.6e l0 ° 6.6e l01 
C" 0.5e-' 2 ' 5 0.6e- l ° 4 1.9 e -i1 - 3 0.6 e - l °' 2 OJe" 10 - 4 1.7e _i1 ' 3 

Pe^ 5.8e' 31 5.8e~' 31 5.4e~' 30 5.8e~' 31 5.8e~' 31 5.4e~' 30 

T 24.3e i0 '° 23.5e l0 ° 23.1e i0 ° 23.6e" iai 23.5e l0 '° 23.2e i0 ° 
P 4.7e~ l0A 6.5e" 10 - 4 6.3c" 10 - 3 6.7 e"' 03 5.7e- l °' 4 5.6e~ l0A 
C 0.8e l26 2.2e'°- 2 4.86- 11 - 1 2.3e l04 2.2e l0 - 2 4.3e- L1 

P e » 0.7e 8 °'° 0.7e l0 -° 0.7 e - !01 0.7e ,;0 '° 0.7e l0 -° O^e" 10 ' 1 



TABLE V: Topological amplitudes in units of 10 5 GeV for the B — > irK, tttv decays in the NDR scheme. 



Data LO LOnlowc +VC +QL +MP +NLO 





0.31 ±0.26 


0.70 


0.73 


0.74 


0.73 


0.73 


„ ? .+0.02 (+0.01) 
-0.03 (-0.01) 




-0.50 ±0.12 


-0.34 


-0.49 


-0.47 


-0.51 


-0.41 


n 49 +1.00 (+0.05) 
U '^ Z -0.56 (-0.05) 



TABLE VI: Mixing-induced CP asymmetries in the NDR scheme. 



results of S„af( s are close to that extracted from the 6 — > ccs decays, S ces — sin(20i) « 0.685, as shown in 
Table IVI1 On the contrary, P is of 0(A) of T in the B° — > tt^tt^ decays, such that a larger deviation of S nn from 
S c5s is found. The LO PQCD results of S vn are consistent with the data, but those of S^o Ks are not. Moreover, 
PQCD predicts AS 7T o Ks = S n o K —^ S rF , K > 0, opposite to the observed values. This result is in agreement with 
those obtained in the literature |l5l I37L l38|. Hence, the measurement of the mixing-induced CP asymmetries 
in the penguin-dominated modes provides an opportunity of discovering new physics. Currently, the data of 
S-k°Ks s t n l suffer significant errors. On the other hand, the NLO corrections and the theoretical uncertainty, 
which concern the allowed range of the PQCD predictions, need to be analyzed. A more clear picture will be 
attained, after we complete these analyses. 



III. NEXT-TO-LEADING-ORDER CORRECTIONS 



We explain the consistent power countings in a s and in large logarithms, before computing the NLO correc- 
tions. A PQCD formula of leading-power in l/m& is written symbolically as 

exp[ 7 (°)(a s )i wc ] <g> exp[r(°)(Q s )L| + T^(a 2 s )L 2 s ] ® exp[ 7 ( )(a s )L flG ] <g> H^(a s ) <g> <£(/z ) , (22) 

where the first, second, and third exponentials represent the Wilson coefficient, the Sudakov factor, and the 
RG factor [8fll |. with the notations Lwc = ln(ffiff/t), Ls = ln(xP6), and Lrg = ln(t&), xP being a fractional 
parton momentum, t ~ \Mmj a characteristic hard scale, and b the conjugate variable to the parton transverse 
momentum kx, and 7, T, and j q the corresponding anomalous dimensions. The RG factor governs the evolution 
from t down to 1/6. The evolution from 1/6 down to the cutoff /iq, which characterizes the meson distribution 
amplitude 0, has been neglected. This formula is complete at LO, since the hard kernel H is evaluated to 0(a s ), 
and at next-to-leading logarithm (NLL), since the Wilson coefficient, the Sudakov factor, and the RG factor 
have been resummed up to the next-to-leading logarithms a s Lwc, ot s Ls, and a s Lnc, respectively (a s Lg is 
the leading logarithm). In all our previous works we used the one-loop running coupling constant a s , which 
is, strictly speaking, a NLO effect. This effect takes into account the potential large NLO corrections to the B 
meson transition form factors through the BLM procedure (see Introduction). 
Next, we add subleading corrections to Eq. Q22[l. which include 

1. H^(a s ) — ► H^°'(a s ) + H^ x \a 2 ). This is what we are going to do in this section, where the NLO hard 
kernel H^ 1 ' contains the vertex corrections, the quark loops, and the magnetic penguin. 



9 



2. exp[7^ (a s )Lwc] — ► exp^ ) (a s )Lwc + J^ 1 '( c ^)^ j wc]- The LO Wilson coefficient is replaced by the 
NLO one, for which the corresponding anomalous dimension is calculated to two loops: 7^°^(a s ) — * 
7 < -°' ) (a s )+7^ 1 ' (a 2 ). According to our counting rules, the NLO anomalous dimension leads to the summation 
of the next-to-next-to-leading logarithm (NNLL) a 2 Lwc- 

3. exp[T^(a s )L 2 s +T^{a 2 )L 2 s ] -> exp[r^(a s )L| + T^\a 2 s )L 2 s + T^{al)L 2 s }. This means the accuracy of 
the summation up to NNLL (a 3 s L 2 s ). Unfortunately, it requires a three- loop evaluation of the corresponding 
anomalous dimension for the Sudakov factor, which is not yet available in the literature. 

4. exjp[y q°\a s ) Lug) — * sxp[yg {a s )LnG + (o 2 )Lhg]. Since Lug an d Lg are of the same order of mag- 
nitude, and the NNLL Sudakov resummation is not available, this NNLL piece of subleading corrections 
(a 2 Lug) can not be included consistently. 

The power countings in a s and in various large logarithms are independent in principle. Based on the above 
classification, we shall extend Eq. (|22|) by considering the subleading corrections from the first and second 
pieces. With the one-loop running a s , the NLO corrections to the hard kernel are complete (assuming that the 
corrections to the form factors have been minimized by our choice of the hard scale) . It is not necessary to adopt 
the two-loop a s as in [33|, whose effect is next-to-next-to-leading-order (NNLO). Because of Lwc ^> Ls,Lrg, 
the NNLL term a 2 L\vc is much more essential than those from the third and fourth pieces. The LO PQCD 
results for the B — > nK, -kit decays from using the NLO Wilson coefficients have been listed in Tables ITlIIVIl 
When investigating the NLO corrections from the vertex corrections, the quark loops, and the magnetic penguin 
to the hard kernel below, we shall always use the NLO Wilson coefficients. After obtaining the decay amplitudes 
A(B — > /) up to NLO, we employ Eq. (fT4"|) to evaluate the corresponding decay rates. 



A. Vertex Corrections 



It has been known that the vertex corrections, reducing the dependence of the Wilson coefficients on the 
renormalization scale /z, play an important role in a NLO analysis. Since the nonfactorizable contributions are 
negligible [40(. we concentrate only on the vertex corrections to the factorizable amplitudes. For charmless B 
meson decays, these corrections do not involve the end-point singularities from vanishing momentum fractions 
in collinear factorization theorem (QCDF ). Therefore, there is no need to employ kx factorization theorem 
(PQCD @> Sim ^1^2]) here. This claim can be justified by recalculating one of the nonfactorizable amplitudes, 
A4 e i, for the B — > ttK decays in collinear factorization theorem, which is also free of the end-point singularity. It 
is found that the results for M e 4 from the two calculations (with and without the parton transverse momentum 
kx in the kaon) differ only by 10%. For more detail, refer to Appendix. After justifying the neglect of the 
parton transverse degrees of freedom, we simply quote the QCDF expressions for the vertex corrections. An 
important remark is that the li ght quark from the b quark transition is assumed to carry the full momentum of 
the associated meson in QCDF |2J|. Strictly speaking, this light quark carries the fractional momentum, whose 
dependence should appear in the PQCD formalism for the vertex corrections. Because it is indeed an energetic 
quark, the assumption is reasonable. 

The vertex corrections modify the Wilson coefficients in Eq. 112fl into [2j| 

An N c 
azip) -> a 2 (/i) + -^j— C F — —V 2 {M) , 

0,0*) - + ^C F ^±^y,(M) , i = 3 - 10, (23) 

An N c 

with M being the meson emitted from the weak vertex. For the B — > irK decays, M denotes the kaon for the 
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vertex functions Vi,4,6,8,io and the pion for 1^2,3,5, 7,9 



In the NDR scheme Vi (M) are given by [25| 




dx 4>m{x) g(x) 



for i 



1-4,9,10, 



Vi(M) 



< 



dx (f>ff(x) g(l - x) , for i 



5,7 



(24) 




v JM Jo 

where /m is the decay constant of the meson M, and (j) M {x) [<j^{x)\ the twist-2 (twist-3) meson distribution 
amplitude given in Sec. IIVI x being the parton momentum fraction. The hard kernels are 

g( X )=3( 1 -^ln X -in) 



+ 2Li 2 (:c)-hr x -\ (3 + 2i 7r) In x — (x «-> 1 — x)\ , (25) 

1—x 

h(x) = 2Li 2 (x) - ln 2 x - (1 + 2i7r)ln x — (x «-» 1 — x) . (26) 



The expressions of Vi (M) in the HV scheme can be found in 43] . The factorization formulas for the various 
B — > irK, tttt decay amplitudes are still the same as in Tables HI and UTI 

The dependence of the Wilson coefficients aj(/x) on the renormalization scale y, modified by the vertex cor- 
rections is exhibited in Fig. ^ for both the real and imaginary parts. It is found that the /j, dependence of most 
of Oi is moderated by the vertex corrections (with the generation of the imaginary parts) . The fj, dependence 
of (16,8 is, however, not altered. It has been known that their dependence will be moderated after combined 
with the fi dependence of the chiral scale ttiok (m) associated with the kaon The most dramatic changes 
arise from 02,3,10- Due to the smallness of 0,3 (aio) compared to the Wilson coefficient 04,6 (09) for the QCD 
(electroweak) penguins, the only significant effect appears in the color-suppressed tree amplitude C , which is 

governed by a 2 . For other aj, the vertex corrections amount only up to 70% at the scale /i ~ ^Amj ~ 1.5 
GeV. The above observation is manifest in Table IVI most of the topological amplitudes for the B — > ttK, tttt 
decays change a little, while C and C are enhanced by factors of 3 and 2 (viewing the values in the columns 
LOnlowc and +VC ), respectively. 

It is then understood that the B — > ttK branching ratios, dominated by the penguin contributions from 04,6, 
vary only slightly under the vertex corrections, as indicated in Table ITlTl However, the direct CP asymmetries 
of the B ± vr ^ and B° -> tt°K° modes, related to C, are modified significantly, as shown in Tables llVl 
Acp(i? ± — * tt^K^) has increased from —0.06 to —0.01, and Acp(B° — > ir°K°) = A^o Ks has decreased from 
0.00 to —0.07. Acp(B° — > tt^K^), determined solely by the color-allowed tree amplitude T' , does not change 
much. The effect from the vertex corrections on the LO PQCD predictions for the B — > tttt decays can also 
be understood by means of the enhanced color-suppressed tree amplitude C: the B° — > tt°tt° branching ratio 
increases by 30%, and the direct CP asymmetry changes from —0.34 to +0.65. The sign flip of the direct 
CP asymmetry is attributed to a huge change of the strong phase of C caused by the vertex corrections. The 
predicted B° — > tt^tt^ and B ± — > tt^tt branching ratios, to which C remains subdominant, decrease only a bit. 
The NLO effect, though increasing \C\ by a factor of 2 , is not enough to resolve the B — > tttt puzzle. Perhaps, 
the penguin amplitude is also larger than expected |3(J,|44|. Nevertheless, the vertex corrections do improve the 
consistency between the theoretical predictions and the experimental data of the B — > tttt decays. 

Though the vertex corrections have been included in QCDF JgBJ, they do not help resolve the B — * ttK puzzle. 
We neglect the electroweak penguin P' ew for convenience in the following explanation. Table shows that the 
penguin amplitude P' is in the second quadrant, and the color-allowed tree amplitude T" is roughly aligned with 
the positive real axis. The color-suppressed tree amplitude C is enhanced by the vertex corrections, and becomes 
almost imaginary. It then orients the sum T' + C into the fourth quadrant, such that T' + C and P' more 
or less line up (and point to the opposite directions) . This is the reason the mag nitude of A CP (B ± -> tt ^), 
proportional to the sine of the angle between T' + C' and P', becomes smaller in PQCD. The situation in QCDF 
is different, where P' is preferred to be in the third quadrant 40] . That is, the predicted Acp{B° — > tt^K^) 
has a wrong sign. Then the modified C, still orienting T' + C into the fourth quadrant, can not reduce the 
magnitude of Aqp{B ± — * tt^K^. The three types of NLO corrections considered here have been extended up 
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FIG. 1: Real parts of a,i for the B — > tt/C decays without the vertex corrections (dotted lines) and with the vertex 
corrections (solid lines) , and imaginary parts with the vertex corrections (dot-dashed lines) in the NDR scheme. 



to 0(alPo) in QCDF recently |4^, which, however, make the QCDF predictions for AcpiB^ — > tt^K^) more 
deviate from the data. Another 0{a 2 s ) piece from the b — > sg*g* transition was included into QCDF ^(|, which 
enhances the B — > 7rii" branching ratios, but leaves their direct CP asymmetries intact. The S — > 7rif puzzle 
can not be resolved in SCET either the leading SCET formalism requires the ratio C'/T' to be real, so 
that C, being parallel to T', can not orient the sum T" + C into the fourth quadrant, and that the magnitude 
of Acp{B ± — > 7r° J ftT ± ) remains large. 

We have found that the color-suppressed tree amplitude C could be enhanced few times by the vertex cor- 
rections in the standard model. It is then worthwhile to investigate whether the mixing-induced CP asymmetry 
S^oks m tnc B — > tt°Ks decays deviates from S c5s substantially under a large C according to Eq. ljTT|) . A 
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similar investigation of the large C effect applies to S^v in the B° — ► -k^tt^ decays according to Eq. (|20l) . The 
results are collected in Table IVII which indicates that the deviation is still small and positive. It is known 
that the leading deviation caused by C is proportional to cos(<5c" — #P')j if neglecting P' ew . Because the vertex 
corrections also rotate the orientation of C", it becomes more orthogonal to P' as shown in Tabic Fvl and AS 7r o Ks 
is not increased much. The mixing- induced CP asymmetry S^, depending only on T and P, is insensitive to 
the vertex corrections, which mainly affect C. 



B. Quark Loops 



For the B — ► TtK and B — > irir decays, the dominant penguin amplitude P' ~ | Vt6 V^* | C4 and tree amplitude 
T ~ \V u bV* d \C2 are both of 0(A 4 ) 0|- Hence, the charm-quark loop amplitude, proportional to a s |Vd,V^,* IC2 ~ 
a s X 2 in the former and to a s |y c 6V c ^|C2 ~ a s X 3 in the latter, could be an important source of NLO corrections. 
Its effect is expected to be larger in the B — > irK decays. On the other hand, the up-quark loop amplitude, 
proportional to a s |V^,V^ s |C2 ~ a s A 5 for B — > ttK, seems to be negligible. For B — ► 7T7t, the up-quark loop 
amplitude, proportional to a^T^V^IC^ ~ a sA 4 J3j might be comparable to the charm-quark one. Therefore, 
we shall include both quark loops in the following analysis. For completeness, we shall also include the quark- 
loop amplitudes from the QCD penguin operators, whose contributions are proportional to a s |Vtf,VJ* |Cj ~ a s A 4 , 
i = 3, 4, 6. They have the order of magnitude the same as or larger than the up-quark one, and can influence the 
direct CP asymmetries of the B — > ttK modes. The quark loops from the electroweak penguins will be neglected 
due to their smallness. Note that the CKM factors of these corrections differ among the loop amplitudes and 
between the b — > s(d) and b — > s(d) transitions. 

For the b — > s transition, the contributions from the various quark loops are given by 



H. 



ctr 



q—u,c,t q' * 



(27) 



I 2 being the invariant mass of the gluon, which attaches the quark loops in Fig. El For the b —> d transition, the 
quark-loop corrections are obtained by substituting d for s in Eq. I|27|) . The functions C^ q \n,l 2 ) are written as 



G^(^l 2 )-^ 



<? 2 (M) 



(28) 



for q = u, c, and 



G^^,l 2 )-± 



C 3 (aO + E G^")(/x,Z 2 )[C 4 (m) + C6(m)] 

q"=u,d,s,c 



(29) 



The constant term —2/3 in the above expressions arises from the Fierz transformation of the four-fermion 
operators in D dimensions with the anti-commuting Dirac matrix 75 in the NDR scheme. The contribution 
proportional to the Wilson coefficient C5, being purely ultraviolet, should be combined with that from the 
magnetic penguin to form the effective Wilson coefficient Cs g + C5 [IJ. Since our characteristic hard scale is 
of order ^Anij, ~ 1.5 GeV, the b quark is not an active one, and does not contribute to Eq. (|29|l . Except this 
difference, our expressions are basically the same as in [2^ • 

The function G^ c \n,l 2 ) for the loop of the massive charm quark is given by 



G (C W 2 ) 



1 2 
mi 



-4 / dx x(l — x) In 
/o 



x(l - x)l 2 



(30) 



m c being the charm quark mass, whose real and imaginary parts are 



2 (h 



3 V 3 



im 2 , mf 
—rrr- - In —k- 
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FIG. 2: Quark-loop amplitudes. 
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i-- 
-i 



4 m ^ 



-2 1 



1-^ln- 



i+i i- 
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ImG<W) = ?(l 



2m^ 
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(31) 



(32) 



respectively. For the loops of the light quarks u, d, and s, we have the expressions similar to Eq. I|3U|) but with 
m c being replaced by m u , rrid, and m s , respectively. Because their contributions are insensitive to the light 
quark masses, we simply adopt the same mass m for the three quark loops. Varying m from m u — 4.5 MeV to 
m s w 100 MeV, the branching ratios change by less than 1%. 

To picture the quark-loop effect, we display in Fig.^the dependence of [VqbV*, d - ) /(VtbV t *, d s)]C t - q \ q = u,c,t, 



on the renormalization scale /x for a given I 



j/4 in the NDR scheme. The real part of the up-quark 



loop is indeed negligible compared to that of the charm-quark loop in the b — > s transition as indicated in 
Fig. Eta)- However, in the other transitions described by Figs. EJb)-|3{d), the up- and charm-loop corrections 
are comparable as argued above. The quark loops from the QCD penguin operators are in fact essential. 
Figures GHa) and|!Jc) [and also Figs.|3{b) and^d)] imply that the weak phases cause different fi dependences 
between the b — > s and b — * d transitions in the cases of the up and charm loops, but not in the case of the 
QCD-penguin loops. 

The quark- loop amplitudes depend on the gluon invariant mass I 2 , which is assumed to be an arbitrary 
constant (I 2 ) in FA. Since the topology displayed in Fig. [21 is the same as the penguin one, its contribution can 
be absorbed into the Wilson coefficients 04,6, 



04,6 (m) -> a 4,6 (/■*) + 



a 



.00 



9tt 



E 

q—u,c,t 



V qb V, 



v tb v t 



5£W c( <r) M2)) 



(33) 



ts(d) 



with the other a% unmodified. The resultant values of 04^ at /x = 1.5 and 4.4 GeV are listed in Table IVTT1 
As fi — 1.5 GeV, the quark-loop corrections do not change a&fi much for b — > s and b — > s, while they are 
destructive (constructive) to 04,6 for b — * d (b —> d). As /x = 4.4 GeV, these corrections are always constructive 
for the different b quark transitions. Besides, the quark-loop corrections are mode-dependent. For example, 
they cancel between the uu and dd components of 7r° = (uu — dd)/y/2 in the i? ± — * 7r ± 7r° decays, but do not 
in others. 

The assumption of a constant gluon invariant mass in FA introduces a large theoretical uncertainty as making 
predictions. In the more sophisticated PQCD approach, the gluon mass is related to the parton momenta 
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FIG. 3: Quark-loop contributions to the b — > s[(a) and (b)] and b — > d[(c) and (d)] transitions for I 2 — m 2 B /A with the 
solid, dotted, and dot-dashed lines corresponding to the up-quark, charm-quark, and QCD-penguin loops, respectively. 





LOnlowc +QL {b -> s) +QL (6 -> d) +MP 


a 4 (1.5 GeV) 
a 6 (1.5 GeV) 


-0.0601 -0.0659 - i 0.0152 -0.0500 - i 0.0131 -0.0492 
-0.0952 -0.1010 - i 0.0152 -0.0850 - i 0.0131 -0.0843 


a 4 (4.4 GeV) 
a 6 (4.4 GeV) 


-0.0336 -0.0545 - i 0.0048 -0.0454 - i 0.0036 -0.0279 
-0.0428 -0.0637- i 0.0048 -0.0546 - i 0.0036 -0.0371 




LOnlowc +QL (6 -> s) +QL (6 -> d) +MP 


a 4 (1.5 GeV) 
a 6 (1.5 GeV) 


-0.0601 -0.0646 - i 0.0150 -0.0804 - i 0.0180 -0.0492 
-0.0952 -0.0997 -i 0.0150 -0.1155 - i 0.0180 -0.0843 


a 4 (4.4 GeV) 
a 6 (4.4 GeV) 


-0.0336 -0.0537- i 0.0047 -0.0628 - i 0.0065 -0.0279 
-0.0428 -0.0630 - i 0.0047 -0.0720 - i 0.0065 -0.0371 



TABLE VII: a 4j 6 including the quark loops and the magnetic penguin for I 2 = m 2 B /4, in the NDR scheme. 



unambiguously (see Appendix). Due to the absence of the end-point singularities associated with I 2 , 1' 2 — > in 
Figs.[2Ja) and|2Jb), respectively, we have dropped the parton transverse momenta kr in I 2 , 1' 2 for simplicity. 
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The amplitudes in Eq. 10 become 



(t) 

irK ' 



a(u,c) JulM At) M) _ jAt) 

.(u,c) A (u,c) ,J_ M (u,c) ,(t) ,{t) l —KA 

Au,c) A (u : c) , M (u,c) At) At) _A / |(*) 

A (u,c) A {u,c) LuM A {t) -> 4 (t) x_yM 

Au,c,t) _ ^(u.c.t) 



(34) 



where M* and .M^ denote the up-, charm-, and QCD-penguin-loop corrections, respectively, and the 

minus sign for the final state tt°K° comes from the (id component in tt°. The factorization formulas for M^^'^ 

and M^ '^ are presented in Appendix. 

As indicated in Eq. (|33[l . the quark-loop corrections affect the penguin contributions, but have a minor impact 
on other topological amplitudes. This observation is clear in Table [V] \P'\ (\P\) has increased from 50.6 to 
52.1 (6.5 to 6.7) in the NDR scheme. Since the B — > ttK decays are penguin-dominated, their branching ratios 
receive an enhancement (see Table ITTTf) . The increase of the branching ratios then reduces the magnitude of 
the direct CP asymmetries in the B — > ttK modes slightly as shown in Table llVl It is also easy to understand 
the insensitivity of the mixing-induced CP asymmetry S^oj^g to the quark- loop corrections (see Table IVT|) . 
viewing the small change in the dominant amplitude P' in Eq. 1)17(1. On the contrary, the penguin contribution 
is subdominant in the B — > irir decays, so the branching ratios do not vary much. However, the direct CP 
asymmetries Acp(B° — > tt^^) and Acp(B a — > 7r°7r°), and the mixing- induced CP asymmetry S^, directly 
related to the penguin amplitude, change sizably. 



C. Magnetic Penguins 

We then discuss the NLO corrections from the magnetic penguin, whose weak effective Hamiltonian contains 
the b — ► sg transition, 

H eS = -^=V tb V* s C 8g 8g , (35) 

with the magnetic penguin operator, 

°s a = ^rn b SiCT^{l + lb )T«G a ^b , (36) 

i, j being the color indices. The Hamiltonian for the b — > d transition is obtained by changing s into d in 
Eq. I|35|l . The topology of the magnetic penguin operator is similar to that of the quark loops. If regarding 
the invariant mass I 2 of the gluon emitted from the operator 8g as a constant (l 2 ), the magnetic-penguin 
contribution to the B — > irK, irir decays can also be included into the Wilson coefficients, similar to Eq. (|33|l . 



a ifi {n) -> 04,8 (/*) - 7fT^ C tg 



0*) > (37) 



with the effective Wilson coefficient C 8g — C 8g + C5 |3jJ. The resultant Wilson coefficients 04,6 (/x) for /x = 1.5 
and 4.4 GeV have been presented in Table IVIII The cancellation between the real parts of the quark-loop 
corrections and of the magnetic penguin is obvious, except in the case of the b — > d transition for ji = 1.5 GeV. 
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In the PQCD approach the gluon invariant mass I 2 is related to the parton momenta, such that the corre- 
sponding factorization formulas involve the convolutions of all three meson distribution amplitudes. Because 
the nonfactorizable contributions are negligible, we calculate only the magnetic-penguin corrections to the the 
factorizable amplitudes, which modify only A* in Eq. (0: 

A (t) -> A {t) -M (9) 

A [t) -> A [t) - —M {9) 

A {t) -> A {t) -M {a) 

TT^TT 7T ' 7T 7T7T ■ 

At) A*) 

4to - 4to - ^JU« ■ (38) 

The explicit expressions for the magnetic-penguin amplitudes M^fjc and mS^I are referred to Appendix. Since 
an end-point singularity arises, as the invariant mass I 2 approaches zero, we have employed the kr factorization 
theorem, i.e., the PQCD approach in this case. 

The effect of the magnetic penguin is just opposite to that of the quark-loop corrections as indicated in 
Tables Hill VI it decreases all the B — > irK, nn branching ratios, except those of the tree-dominated B° — > 7r T 7r ± 
and B + — > ^ir modes, and intends to increase the magnitude of most of the direct CP asymmetries. The 
mixing-induced CP asymmetry S^oks is stable under the magnetic-penguin correction for the same reason. The 
magnitude of S^n decreases due to the smaller penguin pollution. Because the quark-loop corrections are smaller 
than the magnetic penguin, the pattern of their combined effect is similar to that of the latter. In summary, 
the above two pieces of NLO corrections reduce the LO penguin amplitudes by about 10% in the B — ► irK, ttit 
decays, and the B — > ttK and B° — ► ir°ir° branching ratios by about 20%. The direct CP asymmetries are not 
altered very much, which are mainly affected by the vertex corrections, as shown by the similarities between 
the columns +VC and +NLO in Table ITVI 



IV. THEORETICAL UNCERTAINTY 



In this section we explain in detail how to derive the results in Tables ITlIIVIl and discuss their theoretical 
uncertainty. The PQCD predictions depend on the inputs for the nonperturbative parameters, such as decay 
constants, distribution amplitudes, and chiral scales for pseudoscalar mesons. For the B meson, the model wave 
function has been proposed in |21j: 



(j) B {x,b) = N B x 2 (l - x) 2 exp 



1 / xrriB 

2 I 



co 2 B b 2 



(39) 



where the Gaussian form was motivated by the oscillator model in |48|. and the normalization constant Nb is 
related to the decay constant fs through 

dxct)B(x,b = 0) = —^=- (40) 



2y/2N c 

There are certainly other models of the B meson wave function available in the literature (see 0,|H(j). It has 
been confirmed that the model in Eq. I|39|) and the model derived in [5l| with a different functional form lead 
to similar numerical results for the B — > 7r transition form factor 1521. 



The twist- 2 pion (kaon) distribution amplitude <fin(K) 



parameterized as 

,(*) = 



, and the twist-3 ones <l>n(K) an< ^ ^tt(k) nave been 



Itt(K) 



6x(l 



1 + al (K) Cl /2 {2x - 1) + al (K) Cl /2 {2x - 1) 



n(K) r ,3/2 



Cf\2x-l) , (41) 
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\(K) 



(,;) 



Uk) ^ + (30^-^)^(2,-1) 

- 3 {^3 + ^W 1 + ^1 {K) )) Cl /2 (2x - 1) 



(1 - 2x) 



1 + 6 577. 



p£ w - \pI{k)4 k) ) (1 - 10s + 10^) 
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(42) 
(43) 



with a\ = 0, the mass ratio p n (K) — ( m u + m d{s))/ m Tr(K) — m Tr(K)/ m 0Tr(K) an< i the Gegenbauer polynomials 



Cl'\t) = \{U 2 -l) , C'l /2 (t) = i(3-30i 2 



35 1 4 



C? /a (t) = 3t, 



C 2 3/2 (*) = |(5t 2 -l) 



C 4 3/2 W = ^ (1 - 14t 2 + 21i 4 ) 
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(44) 



In the above kaon distribution amplitudes the momentum fraction x is carried by the s quark. For both the 
pion and kaon, we choose 7/3 = 0.015 and L03 = —3 [&|. Because we did not employ the equations of motions for 
the twist-3 meson distribution amplitudes |25| . we are allowed to include the higher Gegenbauer terms, which 
are in fact important. However, we drop the derivative term with respect to the transverse parton momentum 
kx in <jr[(K)- It has been observed that the contribution from this derivative term to the B — ► 7r form factor is 
negligible |54| . 

In our previous works we adopted the models of the pion and kaon distribution amplitudes derived from QCD 
sum rules in :53|. Fixing the B meson decay constant fs ~ 190 MeV from lattice QCD (see H3), the shape 
parameter of the B meson wave function was determined to be ujb ~ 0.4 GeV [2lJ from the B — > tt form factor 
Ff 7r (0) w 0.3 in light -cone sum rules j5|| |57|]- The chiral scales were chosen as w 1.3 GeV for the pion 
and m,QK ~ 1-7 GeV for the kaon 0- The renormalization scale /i was set to the off-shellness of the internal 
particles, consistent with the BLM procedure. The resultant PQCD predictions have been confirmed by the 
observed B — > ttK branching ratios and B° — > ir^K^ direct CP asymmetry. The consistency indicates not only 
that the above inputs are reasonable, but that the short-distance QCD dynamics has been described correctly 
in PQCD. 

In this paper we employ the updated models of the pion and kaon distribution amplitudes in |58| . Since the 
updated Gegenbauer coefficient aj = 0.115 is smaller than the previous one 0.44 for the twist- 2 pion distribution 
amplitude [53, ^^^(0) reduces compared to that obtained in [24|. To compensate this reduction, we increase 
the B meson decay constant up to Jb = 210 MeV, which is consistent with the recent lattice result |59|. in order 
to maintain the B — > itK, 7T7t branching ratios. For the same reason, the penguin annihilation amplitudes, which 
involve the ?t-K or tt-tt time-like form factor, decrease. The magnitude of the resultant direct CP asymmetries 
of the B — > ttK, 7T7T decays, which is not compensated by the overall decay constant fs, then becomes smaller 
than in [24[ as shown in the column LO of Table llVl The smaller B° — > ir^K^ direct CP asymmetry is in 
better agreement with the data, implying that the data could be covered by the theoretical uncertainty at LO 
of PQCD. 

All the above nonperturbative inputs suffer uncertainties, and it is necessary to investigate how these un- 
certainties propagate into the predictions for nonleptonic B meson decays. Here we shall constrain the shape 
parameter ujb and the Gegenbauer coefficients of the twist-2 pion distribution amplitude 0^ using the experi- 
mental error of the semileptonic decay B — > ixlv. The sufficient uncertainties will be assigned to the Gegenbauer 
coefficients of the twist-2 kaon distribution amplitude (f>j^. The other inputs, such as the B meson decay con- 
stant, the twist-3 distribution amplitudes, and the chiral scale associated with the pion and the kaon will be 
fixed. On one hand, the considered sources of theoretical uncertainties have been representative enough. On 
the other hand, it is impossible to constrain all the inputs with the currently available data. 

The spectrum of the semileptonic decay B — > ttIv in the lepton invariant mass q 2 has been measured |60j : 



f*(dr/dq 2 )dq 2 

rtotal 



(0.43 + 0.11) x 10" 



(45) 



with the total decay rate r to tai = (4.29 ± 0.04) x 10 GeV 61]. Assuming that the above error is uniform 
in the region < q 2 < 8 GeV 2 , we derive the uncertainty A of (dr/c?(7 2 )| g 2 =0 by solving the equation 8A = 
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0.11 x 10 4 r to tai, where we take only the central value of r to tai for simplicity. With the allowed range of 
\V u b\ — (3.67 ± 0.47) x 10~ 3 |§l|, A is translated into the uncertainty of the B — ► tt form factor, 

Ff^O) = 0.24 ±0.05 , (46) 

whose central value comes from our choice of the inputs. Equation (|46Jl is consistent with 0.23 ± 0.04 extracted 
in j62 from a global fit to the above CLEO data, lattice QCD results of F^ 7r (q 2 ), etc. A numerical analysis 
indicates that F+ n (0) is more sensitive to lob than to the Gegenbauer coefficients of <fi^. 
Therefore, we propose 

1. the shape parameters for the distribution amplitudes, 

uj b = (0.40 ± 0.04) GeV , a£ = 0.115 ± 0.115 , a\ = -0.015 , 

af = 0.17 ±0.17 , af = 0.115 ±0.115, of = -0.015 , (47) 

that is, the Gegenbauer coefficients can vary by 100%. We do not consider the uncertainty from the 
coefficients a\ and af , to which our predictions are insensitive. Note that the first Gegenbauer coefficients 
af ps 0.10 ± 0.12 and af » 0.05 ± 0.02 have been found to be smaller in [&| and respectively. A 
hint on the effect from the evolution of the meson distribution amplitudes from 1/b down to the cutoff /i 
(see Sec. Ill) can also be obtained through the above variation of the Gegenbauer coefficients. 

2. the CKM matrix elements, 

V ud = 0.9734 , V us = 0.2200 , \V uh \ = (3.67 ± 0.47) x 10~ 3 , 

V cd = -0.224 , V cs = 0.996 , V cb = 0.0413 , (48) 

where we consider only the representative source of theoretical uncertainties from \V u b\ |6lj |. This source 
is essential for estimating the uncertainty of the predicted direct CP asymmetries. V c b — (41.3 ± 1.5) x 
10 -3 JH3 has a smaller uncertainty, and the other matrix elements have been known more precisely. 
The unitarity condition VtbV^,^ = ~VubVus(£\ ~ ^ cb ^*s(d) 1& then employed as evaluating the penguin 
contributions. 

3. the weak phases, 

0! = 21.6° , fa = (70 ±30)° , (49) 

where the range of the well-measured <j>\ with sin(20i ) = 0.685 ± 0.032 j6^ has been neglected, and the 
range of </>3 is hinted by the determinations |65l |66J , 

03 = 68t" ± 13 ± 11 (Belle, Dalitz) , 
70 ± 31±J§±^ (BaBar, Dalitz) , 
63±il (CKMfitter) , 

64 ± 18 (UTfit) . (50) 

We fix the other parameters, such as the meson decay constants fs = 210 MeV, fx = 160 MeV, f n = 130 
MeV, the meson masses tub = 5.28 GeV, tok = 0.49 GeV, m, = 0.14 GeV, the charm quark mass m c — 1.5 
GeV, and the B meson lifetimes t b « = 1.528 x 10 -12 sec, t b ± = 1.643 x 10~ 12 sec Q. We also fix the chiral 
scales m 07r = 1.3 GeV and m 0K — 1.7 GeV, where the value of m 07r (m 0K ) is close to that (larger than 1.25±0.15 
GeV) obtained in the recent sum-rule analysis The resultant B — > tt, K transition form factors, 

*(0) = 0.24^ , Ff^(0) = 0.36t°;°? , (51) 

respect Eq. (|46|) from the measurement, and are consistent with the estimation from light-cone sum rules |64|. 
If further including the variation of itiok as a source of theoretical uncertainties, we just enlarge the range of the 
B — * ttK branching ratios, but not of the other quantities. We have tested the dependence of our predictions 
on the cutoff /io, which is found to be weak. 
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FIG. 4: Dependence of R, R c , and R„ on 4>z from NLO PQCD with the bands representing the theoretical uncertainty. 
The two dashed lines represent la bounds from the data. 



The above inputs lead to Tables ITlIIVIl where the theoretical uncertainties are displayed only in the columns 
+NLO. The errors (not) in the parentheses represent those from (all sources) the first source of uncertainties. 
It indicates that the nonperturbative inputs, i.e., the first source, contribute to the theoretical uncertainties 
more dominantly in the B — > ttK decays than in the B — > 7T7t decays, because the former depend on the less 
controllable parameters associated with the kaon. We also observe that Acp(B° — > n^K^) and AcpiB^ — > 
it® K~^~ ) always increase or decrease simultaneously, when varying the nonperturbative inputs. Hence, the B — > 
ttK puzzle can not be resolved by tuning these parameters. After including the uncertainties, the predicted 
B° — > 7r 7r° branching ratio and mixing-induced CP asymmetry S 7r o Ks are still far from the data. 

A more transparent comparison between the predictions and the data is made by considering the ratios of 
the branching ratios. The following three ratios of the B — > ttK branching ratios have been widely studied in 
the literature, 

B(B° -> ttFK±) t b+ 

R = lTm± ±WoS = 0-85 ±0.06, 

B\B ± — ► n^K") t b o 

B(B ± -» TT ^) 

** = l l%°z n 7J) =°- 82±0 - 08 > (52) 

whose values are quoted from [lj. We have confirmed that these ratios depend on the nonperturbative inputs 
weakly. Therefore, their deviation from the standard-model predictions could reveal a signal of new physics, 
such as a large electroweak penguin amplitude. Table ITTTl shows that for 03 = 70°, the ratio R increases slightly 
from 0.90 to 0.92, when the NLO Wilson coefficients are adopted, beyond which the various NLO corrections 
do not change R much. The ratio R c (Rn) decreases from 1.20 (1-25) to 1.14 (1.14), when the NLO Wilson 
coefficients are adopted, and settles down at this value as indicated by the column +NLO. The different types 
of NLO corrections cause only small fluctuations. Comparing the columns LO and +NLO, the consistency 
between the PQCD predictions and the data has been improved. 

Varying the weak phase <p3 and the inputs, we find that the PQCD predictions for R and R n are in good 
agreement with the data in Eq. 152|l . which is obvious from Fig. 21 However, the predictions for R n exhibit 
a tendency of overshooting the data, which is attributed to the smaller PQCD results for the B° — > ir°K a 
branching ratio. A smaller Gegenbauer coefficient of 0^ enhances R n . That is, when using the updated pion 
distribution amplitudes from |58j| . the consistency of the predictions for R n with the data deteriorates. A smaller 
u>b enhances R n . This is the reason we do not lower ub in order to compensate the reduction from the smaller 
aj- Note that tuok has an effect on the electroweak penguin amplitude, i.e., on the B° — > tt°K° branching 
ratio. Hence, we have also studied the dependence of R n on the chiral scale vtiqk- A smaller indeed reduces 
R n , but does not help much: choosing tuok =1.3 GeV causes only few percent reduction of R n . It has been 
known that the B° — > ir°K° branching ratio can be significantly increased by rotating the electroweak penguin 
amplitude P' ew away from the penguin amplitude P' (their values in Table IVl are roughly parallel to each other). 
Therefore, we can not rule out the possibility that P' ew acquires an additional phase from new physics effects 
0,|H,|6^|. However, our theoretical uncertainty is representative, and the actual uncertainty could be larger, 
such that the discrepancy is not serious at this moment. We do not discuss the ratios relevant to the B — > mr 
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decays, because the PQCD predictions for the B° — ► 7r°7r° branching ratio are currently far below the measured 
values. 

V. CONCLUSION 

The LO PQCD has correctly predicted the direct CP asymmetry Acp(B° — > ir^K^), but failed to explain 
another one Acp{B ± — > ifiK^) Q- Phenomenologically, the substantial difference between Acp{B ± — > tt^K^) 
and Acp(B° — ► 7r =F A" ± ) has led to the conjecture of new physics @,I3- However, the difference can also be 
attributed to a large color-suppressed tree amplitude C as pointed out in fl3l | . Theoretically, an examination 
of NLO effects is always demanded for a systematic approach like PQCD. Since C itself is a subdominant 
contribution, it is easily affected by subleading corrections. Hence, before claiming a new physics signal in 
the B — > ttK data, one should at least check whether the NLO effects could enhance C sufficiently. This is 
one of our motivations to perform the NLO calculation in PQCD for the B — > ttK, tttt decays here. Another 
motivation comes from the mixing-induced CP asymmetries in the penguin-dominated modes, some of which 
also depend on the color-suppressed tree amplitudes. To estimate the deviation of S n o Ks from S cSs within the 
standard model, one must compute C reliably. 

In this paper we have calculated the NLO corrections to the B — > ttK, tttt decays from the vertex corrections, 
the quark loops, and the magnetic penguin in the PQCD approach. The results for the branching ratios and 
CP asymmetries in the NDR scheme have been presented in Tables IIIIIVII and discussed in Sec. III. It has 
been shown that the corrections from the quark loops and from the magnetic penguin come with opposite signs, 
and sum to about 10% of the LO penguin amplitudes. Their effect is to reduce the B — > ttK branching ratios, 
to which the penguin contribution is dominant, by about 20%. They have a minor influence on the B — > tttt 
branching ratios, and CP asymmetries. The vertex corrections play an important role in modifying direct CP 
asymmetries, especially those of the i? ± — * 7r AT ± , B° — ► tt°K°, and B° — ► 7r°7r° modes, by increasing the color- 
suppressed tree amplitudes few times. The larger color-suppressed tree amplitude leads to nearly vanishing 
Acp(B^ — > 7r°A ± ), resolving the B — » ttK puzzle within the standard model. Our analysis has also confirmed 
that the NLO corrections are under control in PQCD. 

The NLO corrections, though increasing the color-suppressed tree amplitudes significantly, are not enough to 
enhance the B° — > tt°tt° branching ratio to the measured value. A much larger amplitude ratio \C/T\ ~ 0.8 must 
be obtained in order to resolve this puzzle [T^- Nevertheless, the NLO corrections do improve the consistency 
of our predictions with the data: the predicted B° — > 7r ± 7r =F (B° — > 7r°7r°) branching ratio decreases (increases). 

Viewing the consistency of the PQCD predictions with the tiny measured B° — » K°K and B° — > p°p° branching 
ratios, we think that our NLO results for the B — > tttt decays are reasonable. In SCET [20], the large \C/T\ 
comes from a fit to the data, instead of from an explicit evaluation of the amplitudes. The amplitude C was 
indeed found to be increased in SCET by the NLO jet function (the short-distance coefficient from matching 
SCETi to SCETn) [z3, if the parameter set "S4" in QCDF was adopted. The large measured B° -> 7r°7r° 
branching ratio was then explained. However, we emphasize again that the same analysis should be applied to 
the B — ► pp decays for a check. Hence, the B — > tttt puzzle still requires more investigation. 

The tendency of overshooting the observed ratio R n has implied a possible new-physics phase associated 
with the electroweak penguin amplitude P' ew . This additional phase can render P' ew orthogonal to the penguin 
amplitude, and enhance the B a — > 7r° K° branching ratio. We have also computed the deviation AS n o Ks of the 
mixing-induced CP asymmetry, and found that the NLO effects push it toward the even larger positive value. 
Therefore, it is difficult to understand the observed negative deviation without physics beyond the standard 
model. 
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APPENDIX: FACTORIZATION FORMULAS 



We first define the kinematics for the B — » M2-M3 decay, where M 2 (M3) denotes the light pseudo-scalar 
meson involved in the B meson transition (emitted from the weak vertex). In the rest frame of the B meson, 
the B (M 2 , M 3 ) meson momentum Pi (P 2 , P 3 ), and the corresponding spectator quark momentum ki (k 2 , ^3) 
are taken, in the light-cone coordinates, as 

Pi = ^|(1,1,0 T ) , h = (0,XiPf,k 1T ) , 

P 2 = ^(1,0, Or) , k 2 = (x 2 P+,0,k 2T ) , 

P 3 = ^(0,1, Or) , fc 3 = (0,x 3 P 3 ,k 3T ) , (A.l) 

where the light meson masses have been neglected. We also define the ratio r 2 = mo 2 /mB (r 3 = m^/mB) 
associated with the meson M 2 (M 3 ), too 2 ("103) being the chiral scale. 

The factorization formulas for the B — ► M2M3 decay amplitudes appearing in Tables HI and ITT1 arc collected 
below: 

1*1 />oo 

F e 4(a) = 16TrCFm B / dx\dx 2 I bidbib 2 db 2 <j)B{xi,b\) 
Jo Jo 

x I [(1 + x 2 )^fa) + r 2 (l - 2a; 2 ) (tf£ fa) - (fifa))] E e (t)h e (A, B, h,b 2 ,x 2 ) 

+2r 2 ^{xi)E e {t')h e {A\B\b 2 M,x 1 )\ , (A.2) 



Mei(a') 



i>l poo 

Pe6(°) — 32irCpm B / dx\dx 2 / b±db\ b 2 db 2 4>b(xi, &i) 
Jo jo 

x|r 3 [<t>$fa)+r 2 x 2 (^fa)+4^fa))+2r 2 ^fa)] E e (t)h e (A, B,b 1 ,b 2 ,x 2 ) 

+2r 2 r 3 ^fa)E e (t')h e {A',B',b 2 ,b 1 ,x 1 )^ , (A.3) 

/>1 />oo 

P Q 4(a) = lGirC F m B / dx 2 dx 3 / b 2 db 2 b 3 db 3 
Jo Jo 

x j [S3^(X2)^(23) + 2r 2 r 3 2 "(S 2 -) {(^(03) + 0s (*s)) + x 3 (0a (23) " 0a (2a))}] 
xP Q (t)/i e (A,P,6 2 ,6 3 ,a;3) 

- [(1 - x^fa^ffa) - 2r 2 r 3 {-20^ (xj) + x 2 (0f (ica) + ^"(xj))} 0f fa)] 



xE a (t')h e (A',B\b 3 ,b 2 ,x 2 )^ , (A.4) 

/•l />oo 

Pa6(a) = 327rCfTOB / dx 2 dx 3 / b 2 db 2 b 3 db 3 
Jo Jo 

x j [2r 2 ^ fa) (JB3) + rgxg 0^ (35) (0f (Sg) + $ fa))] E a (t)h e (A, B,b 2 ,b 3 ,x 3 ) 
+ [r 2 (l - i 2 ) (0 2 P (^) - 02 3 4 (^) + 2r 3 2 4 (^)0f (a*)] 

xP a (t')/i e (A',P', 63,62,^2)! , (A.5) 

\/2Nc f 1 f°° 
32-kCf m c m 2 B I dx\dx 2 dx 3 I b\dbi b 3 db 3 4>b(xi, h)(f> 3 fa) 

x( [(1- 13)^(35) -r 2 a;2 (0£(35) + <Pxfa))] E'Jt)h n (A, B,b u b 3 ) 
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A B A[ B' 

F e 4,e6 t/x^ULB X\X 2 m,B ^/ximB ^/X\X 2 m,B 

Me4,s5,e6 iy/x2 (x~3 — X\)mB ^JX\X 2 m,B \J X 2 [x\ — ~Xz)mB X\X 2 m,B 

F a 4,a6 i^fxima i^/x^XimB iyfWirriB i^/x^x^mB 

M 0.4,0,5,0,6 

I y/(x! - x 3 )x 2 m B iy/x 2 x i m B yjl - x 2 (x 3 - x±)m B Wx 2 x 3 m B 



TABLE VIII: The invariant masses A, B, A' , and B' in the hard kernels. 



M e6 (a') = 32ttC f 



Ma4(a') 



+ [~(X2 + X 3 )^{x^) + r 2 X 2 (X2) ~ <fi(X2))] E' e (t')h n (A' , B' , bi, b 3 ) 
2N7 



(A.6) 



Nr. 



: 2 
-™>B r 3 



dx\dx 2 dx 3 I b\dbi b 3 db 3 4>b{xx, b±) 



X | [(1 - X 3 )^{- 2 ) (0 3 P (^) - <g(xs)) 

+ra(i - x 3 ) (^m + (^m - 

+r 2 x 2 (^(Sa) - ^(Sa)) (<f>i(m) + 4>I(^))] E' e {t)h n {A,BM,h) 

~ [Xz44&2) (4>i(x3) + 03 (2s)) + r 2^3 (X2) + (xi)) ((& (&) + 03 (^s)) 

+r 2 x 2 {^(M) ~ ^(M)) (0f (Si) ~ 4%(^))] E' e (t')h n (A',B',h,b 3 ) 



(A.7) 



32irC F - 



/2N, 



c ™2 



dx\dx 2 dx 3 I b\dbib 3 db 3 4>B{x\,bx) 



M a6 (a') = 32nC F 



N c 

x { [(1 - x 2 )4>i{xl)^{x^) + r 2 r 3 {(1 - x 2 ) (0 2 p (^) + T (^)) (xs) - <t%&S)) 

+x 3 (^(xl) - (cj>£(xt) +$(X3))}] E' a (t)h n (A,B,b 3 ,h) 

- [x 3 (f>£(x2-)<t>i{x5) + r 2 r 3 |40 P (a^)0f(^) 

-(1 - 13) (^(afc) + 4> T 2 {xl)) (03 P (xi) - <j>T{xi)) 

-X 2 (^(Xl) - </%<&)) (</) P (xi) + ^(5g)) }] ££(i')M^'> S', 63, 61) 

2A7 



(A.8) 



c ™2 



dx\dx 2 dx 3 I bidbib 3 db 3 (f>B(xi,bi) 
Jo 

- ra(l - s 2 ) (0 P (a^) + ^(si)) ^(13) 

+^30 2 4 (5 : 2-) (0 P (x?) - ^(ss)) ] ££(*)M4 S, 63, 61) 
r 2 (l + x 2 ) (^(xj) + ^(3a)) 03 <M) 
+r 3 (-2 + x z )4i{xi) {^{xt} - 4>l{xt)) 1 E' a (t')h n (A', B' , b 3 ,h) \ , 



(A.9) 



where we have adopted the notations ~x 2 = I — x 2 and S3 = 1 — x 3 , and ignored the mass difference between tub 
and m,b. F e Ki and .Me/fi are obtained by choosing M 2 (M 3 ) to be the kaon (pion) in F e i and M. e i, respectively. 

The invariant masses A, B, A' , and B' of the virtual quarks and gluons involved in the above hard kernels 
are functions of Xi, x 2 and x 3 , as in Table IVlIII The hard scales are chosen as 



t = max(v^2|, VF1, 1/60 . 



(A.10) 
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with the index i = 1, 2 for -F e 4,e6, i — 2, 3 for -F a 4,a6i an d i = 1,3 for the nonfactorizable amplitudes. 
The evolution factors Ejp and Ea^ are given by 

E e {t) = a a (i) (i) exp[-5 B (t) - S a (*)] - 

= a s (t)a(t)exp[-S 2 {t)- S 3 {t)} , 

E' e (t) = a s (i) o'(i) exp[-S fl (t) - S a (t) - S 3 {t)]\ b2=bl , 

££(t) = a.(i) a'(t)exp[-S B (t)-S 2 (t)- S 3 (t)]\ ba=ba , 



(A.ll) 



where aS'< represent the combination of the Wilson coefficients appearing in Tables [H and [HI The Sudakov 
exponents associated with the various mesons are written as 



Ssit) = exp 
52 (t) — exp 



-s(x 1 P+,b 1 ) - ^ / -^7 5 (a B (/i)) 

l/bi ^ 



-s(a; 2 P 2 + , 6a) - s((l - x 2 )P 2 + ,b 2 ) - 2 / -^Mf*)) 

i/6 2 M 



(A.12) 
(A.13) 



with the quark anomalous dimension 7^ = —a s /ir. The formula for the exponential S3 is the same as S 2 but 
with the kinematic variables of meson 2 being replaced by those of meson 3. The explicit expression of the 
exponent s can be found in |20l l7lL 172^ . The variable b\, conjugate to the parton transverse momentum kix, 
represents the transverse extent of the B meson. The transverse extents b 2 and 63 have the similar meaning for 
mesons 2 and 3, respectively. For the running coupling constant a s (p), we employ the one- loop expression, and 

the QCD scale Aqq D = 0.250 GeV. The Sudakov exponential decreases fast in the large b region, such that 
the long-distance contribution to the decay amplitude is suppressed. 
The hard functions are written as 



h e (A,B,b 1 ,b 2 ,x i ) = [e(b 1 -b 2 )K (Ab 1 )I (Ab 2 ) 

+6 (6 2 - 61) K (Ab 2 ) I {Abx)} K (Bh) S t ( Xi ) 
h n (A,B,h,h) = K {Ab 3 )[e(bi-b3)K (Bb 1 )I (Bb 3 ) 

+6 (63 - 61) K (Bb 3 ) I (Bh)) , 



(A.14) 



(A.15) 



where St resums the threshold logarithm In x appearing the hard kernels to all orders. It has been parameterized 
as 731 



v^r(i + c ) 



(A.16) 



with c = 0.3. 

The factorization formulas for q = u,c,t, involve the convolutions of all three meson distribution 

amplitudes: 



M 



(?) 

ttK 



n 2 

-16m? 



'2N C jo 



dx\dx 2 dx 3 I bidbib 2 db 2 (f)B(xi,bi) 



x { [(1 + i 2 )^(^)^(Ss) + r w (l - 2x 2 ) (tf (X2) - <PZ(X2)) 0&(35) 
+2rjr#*(5a)0£(5s) + 2r*r K ((2 + x 2 )<f>%(x5) + a&^Sa)) ^(^)] 
xE^(t g ,l 2 )h e (A,B,b 1 ,b2,X2) 

xE^(t' q ,l' 2 )h e (A , ,B , ,b 2 ,b 1 ,x 1 )} , 



(A.17) 



with the evolution factor, 



£ (5) (M 2 ) = [a s (t)] 2 C^(t,l 2 )e^>[-S B (t)-S n (t)} 



(A.18) 
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The hard scales are chosen as 

t q = max ( y/\A 2 \, v^l-B 2 !, \/x^x 3 m B , l/h 



t' q = max(^2|, vl^T, Vks-^lms, 1/6,-) , (A.19) 



with the index i = 1,2. The additional scales \Jx 2 x 3 m B and ^/l 3 ^ — x i \ m B, compared to those appearing in 
Eq. (|A.10I) . come from the gluon invariant masses I 2 = (1 — x 2 )x 3 m 2 B and V 2 — (x 3 — xjjm^ in Figs. |2|[a) and 

ffib). respectively. The formulas for M^ '^ are derived from Eq. (|A.17|) by substituting the mass ratio r n for 
rx, and the distribution amplitudes <f>^ and 4>^ for (f>^ and 0^-, respectively. 
The magnetic-penguin amplitude is written, for the B — > irK modes, as |74j 



M-fc = 16m%— =|= / dxidx 2 dx 3 / foirffoi b 2 db 2 b 3 db 3 4>b[xi, h) 
V 2A/ C Jo Jo 

x { [-(1 - x 2 ) {2^(x^) +r w (3^(^) - ^(a*)) 
4-^X2 (^(x^) 4- 0^(x^))} 0k (xi) 

- r x (1 + x 2 )x 3 ^ (aa) (3<£ K («s) + 0k (*3 ) ) 

- r K r w (l - x 2 ) (0^(xi) 4- <PZ{x5)) (3(/>£(xl) - ^(xj)) 

- r K r 7r x 3 (l - 2x 2 ) (^(x?) - ^(xj)) (3^(xi) + ^(xj))] 
x£; g (< g )ft, g (A, 5, C, 6i, 6 2 , &3, a^) 

- [4r 7r 0^(x 2 ")0j<r(xi') + 2r K r n x 3 (f>% (x^) (30^(xi) 4- 0^(xi))] 

xS g (t;)^(A',S , ,C , ,6 2 ,6 1 ,6 3 ,x 1 )} , (A.20) 

with the evolution factor E g (t), 

E g {t) = [a s (t)] 2 Cf g (t)eMS B {t) - S K (t) - S«(t)] . (A.21) 

Since the terms proportional to rx TV develop the end-point singularities as the invariant mass of the gluon from 
Os g vanishes (x 3 — > 0), we have kept the transverse momentum k 3 x- This is the reason the Sudakov factor 
associated with the kaon appears. The hard function is given by 

h g (A 1 B 1 C, 6i,6 2 ,6 3 ,Xi) 

= -St{x i )K {Bb 1 )K {Cb 3 ) [ ' d6 tan J {Ab x tan 6) J (Ab 2 tan 9) J (Ab 3 tan 9) , (A.22) 

Jo 

with the virtual quark and gluon invariant masses, 



A = ^fx^m B , B = B' = ^x\x 2 m B , C = i\J (1 - x 2 )x 3 m B , 

A' = ^/x\m B , C' — \fx\ — x 3 m B . (A. 23) 

The hard scales t q are the same as in Eq. (IA.19|I with the index i = 1,2,3. 

At last, we present the factorization formula for the nonfactorizable amplitude M. e i with the parton transverse 
degrees of freedom in the kaon being neglected. This formula is employed to justify the approximate equality of 
an amplitude without the end-point singularity evaluated in collinear and kx factorization theorems. Dropping 
the parton transverse momentum k 3 x, the corresponding chansreis to remove the factor m 2 B in Eq. (|A. 6|) and 
the integration J b 3 db 3 . We also replace the hard functions by [75| 

m | / K (AW&x) - K (BW&x) for A^ 2 > \ 
K{A [ '\B^ MM) _> B(l)2 _ A(/)2 ^ ifH w (ytfpr^ K (B(0 bl ) for A^ 2 < ) ' (A ' 24) 
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where and fiW have been defined in Table 171111 Without k 3T , the conjugate variable 62 in the Sudakov 
exponent S v (t) is set to b\. 
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